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Distinct Hagedorn temperatures from particle spectra: 
a higher one for mesons, a lower one for baryonsB 
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We analyze experimental particle spectra and show that the Hagedorn tem- 
perature is significantly larger for mesons than for baryons. The effect can be 
' explained within dual string models: excitations of three strings in the baryon 

produce "faster" combinatorics than a single string in the meson, hence lead to 
a more rapid growth of baryons than mesons. Predictions of other approaches 
for the gross features of particle spectra are also discussed. 
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1. Introduction 

The famous Hagedorn hypothesis Jl|, [|, dating back to pre-chromodynamic 
times of the sixties, states that at asymptotically large masses, m, the density of 
hadronic resonance states, p(m), grows exponentially: 

p{m) ~ exp (Jjpj (1) 

The Hagedorn temperature, Th, is a scale controlling the exponential growth of 
the spectrum. Although the Hagedorn hypothesis has sound thermodynamical 
consequences (one cannot heat-up a hadronic system above this temperature), 
Th should not be immediately associated with thermodynamics. In this talk we 
are concerned with the spectrum of particles per se, as read off form the Particle 
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Fig. 1. Cumulants of meson and baryon spectra, and the Hagedorn-like fit with 
Eqs. (UH), plotted as functions of mass. 

Data Tables ||. In this context the "temperature" Tjj is just a parameter in 
Eq. (0). 

Ever since hypothesis (jl]) was posed, it has been believed that there is one 
universal Hagedorn temperature for all hadrons. Presently available experimental 
data show that this is not the case, as has been pointed out by W. Florkowski 
and WB in Refs. g §. 

This talk has two parts: experimental and theoretical. In the experimental 
part (Sec. |2|) we show how well the Hagedorn hypothesis works even for very 
low masses, and point out the key observation that the mesonic temperature is 
significantly larger from the baryonic temperature. In the theoretical part (Sec. 

we argue that the only framework (known to us) which is capable of producing 
the observed behavior in a natural way are the Dual String Models fTj. In Sec. |] 
we discuss other approaches and more speculative ideas. 



2. Experiment 

2.1. Experimental spectra of mesons and baryons 

In Fig. |l] we compare the cumulants of the spectrum |4j], defined as the number 
of states with mass lower than m. The experimental curve is 

A^xp(m) = ^#;6(m - mi), (2) 
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Fig. 2. Various Hagedorn-like fits, made according to formulas of Table [I]. 



where gi = (2Jj + l)(2/j + 1) is the spin-isospin degeneracy of the ith state, and 
rrii is its mass. The theoretical curve corresponds to 



NtheoAm) 



Ptheor (m')dm', 



where 



PtheorM = /(m)exp(m/T), 



(3) 
(4) 

used in the 
(5) 



with f(m) denoting a slowly- varying function. A typical choice |I] 
plot of Fig. |l], is 

f(m) = A/(m 2 + (500MeV) 2 ) 5/4 . 

Parameters Th and A are obtained with the least-square fit to logA^h e0 r, made 
over the range up to m = 1.8GeV, and skipping the lightest particle in the set. 
Other choices of f(m) give fits of similar quality (see Fig. 0). A striking feature 
of Fig. [I] is the linearity of logiV starting at very low m, and extending till 
m ~ 1.8GeV. Clearly, this shows that (|l|) is valid in the range of available data.El 
However, the slopes in Fig. [I] are different for mesons and baryons. For the 
assumed f{m) of Eq. (||) we get 
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195MeV, T> 



baryon 



141MeV. 



(6) 



This means that T meson > T baryon , and the inequality is substantial! Although 
it has been known to researchers in the field of hadron spectroscopy that the 



Above 1.8GeV the data seems to be sparse and we should wait for this region to be explored 
by future experiments. 
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0.015 
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0.014 
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Table 1. Various Hagedorn-like fits. Rows 1-4 use formulas of Ref. citeSBMl, row 
5 uses a simple exponent, and row 6 uses the scalar string model of Ref. [jTcfl . The 
last two column display the mean suared deviation for the meson and baryon case, 
respectively. 



baryons multiply more rapidly than mesons |§, to our knowledge this fact has 
not been presented as vividly as in Fig. [I]. To emphasize the strength of the 
effect we note that in order to make the meson line parallel to the baryon line, 
we would have to aggregate ~ 500 additional meson states up to m = 1.8MeV as 
compared to the present number of ~ 400. 



2.2. Are we asymptotic? 

An important question is whether the presently available range of masses is 
asymptotic in view of Eq. (H). The answer is no! This is how we can look at this 
question quantitatively. Consider the generic form of the spectrum of Eq. (|j). 
We can rewrite it as 

f(m)e m/T = e lo s/( m )Wr ~ e io S [f(m)+f'(m)Am]+(m+Am)/T _ 



( i i /.M 



Am Am 



const eV I(m> ) = const e T c« . 

where m = fn + Am, and in the range of data m ~ lGeV. We have defined T e g 
as the effective Hagedorn temperature in the (non-asymptotic) region around m. 
The value of T eff follows directly from the data. We have, according to Eq. (7), 

1_ = J_ _ f(m) _ 

T T eff . f(m) ' 

The following statements are obvious: 

• since f'(m) < 0, T < T eff , 

• only at m — > oo we have T = T e g. In the region of data we find significant 
differences between T and T e g. 
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Here is a numerical example. Consider 



= — 2W4> ( 8 ) 



which leads to 

115m 



+ o7=^— (9) 



T T eff . 2 (m 2 + m 2 ) 
Now we take mo = 0.5GeV and m = lGeV and find 

for mesons: T cfr = 311MeV, T = 192MeV (exact fit: 195MeV) 
for baryons: T cS = 186MeV, T = 136MeV (exact fit: 141MeV) 

We conclude that only in the asymptotic region, m >> mo, the choice of 
f(m) is not important. In the region of presently-available data /(m) matters 
very much for the extracted values of the Hagedorn temperature. This simply 
means that we need a theory in order to make quantitative statements! 

The numerical parameters obtained from various choices of the function /(m) 
are collected in Table |l| Figure ^| shows the fits corresponding to the rows 1, 4, 
5 and 6 of Table [T]. Note the fits are very close to each other and the theoretical 
curves are virtually indistinguishable in the region of data. In view of the above 
discussion it makes little sense to treat the Hagedorn temperature as an absolute 
parameter and to quote its value without specifying the model that yields the 
function f(m). 



2.3. Flavor universality 



In Fig. ^ we show the cumulants of particle spectra of a given value of 
strangeness. We can clearly see that the slopes in the figure do not depend 
on strangeness. The meson plot includes various Hagedorn fits of Fig. |2]. The 
two sets of lines are displaced in the m variable by roughly 150MeV, which is the 
difference of the masses on the strange and non-strange quarks. The conclusion 
here is that the addition of the strange quark mass has no effect on the rate of 
growth of the number of states with m. Certainly, we are rediscovering the SU(3) 
flavor symmetry here! 



2.4- Plot in the exponential variable 



We end the experimental part of this talk by showing the same information 
as in Fig. |l], but instead of using logarithmic units on the vertical axes, we 
take exponential units on the horizontal axis. More precisely, we take the fit 
to the spectrum with of the form with the simple exponent (row 5 in Table 
[]]), which leads to the cumulant N(m) = AT(exp(m/T) — 1), where the values 
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Fig. 3. Strange vs. non-strange mesons (a), and baryons (b). 

of A and T result from the least-square fit. Next, we define the variable y = 
AT(exp(m/T) — 1) and plot the cumulants as functions of y. Note that the A 
and T parameters are different for mesons and baryons. Again, the linearity 
of data in the figure is striking. It starts at basically m = 0, and extends to 
m ~ 1.8GeV. The advantage of the plot in Fig. |] to that of Fig. |Tj is that now 
the steps in the experimental cumulant are of a similar size independently of m. 

We conclude this section by stating that the exponential growth of hadronic 
spectra in the region of m up to about 1.8GeV, with T mes > Tb ar , is an experi- 
mental fact. 
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Fig. 4. Cumulants of the meson and baryon spectra plotted in exponential variables. 



3. Theory 

We are faced with two basic theoretical questions: 

1. Why is the spectrum of resonances exponential? 

2. Why do mesons and baryons behave so differently? 

Concerning the first question, let us stress that it is not at all easy to get 
an exponentially rising spectrum of resonances. Take the simplistic harmonic- 
oscillator model, whose density of states grows as m" 1 " 1 , with d denoting the 
number of dimensions. For mesons there is one relative coordinate, hence p ~ m 2 , 
whereas the two relative coordinates in the baryon give p ~ m 5 . Weaker-growing 
potentials lead to a faster growth of the number of states, but fall short of the 
behavior ([!]). We know of three approaches yielding behavior ([!]), both involving 
combinatorics of infinitely-many degrees of freedom. These are the Statistical 
Bootstrap Model [@, 0,0,0, Bag Models |L2|, 0, 0, and Dual String Models 
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The first two, however, lead to the same rate of growth for the mesons and 



baryons. Statistical Bootstrap Models are discussed in Sec. pTT| . In Bag Models 



[12], [13], [14]] the exponential growth of the spectrum is associated with the melting 
out of the vacuum around the bag when the hadron is being excited. Since the 
scales in the Bag Model are practically the same for the meson and the baryon 
(the size scales as the number of constituents to the power 1/4), the Bag Models 
are not capable of answering question 2. On the other hand, the Dual String 
Models is offer a natural explanation of questions 1 and 2. This has already 
been pointed out in Ref. f|. 



3.1. Statistical Bootstrap Models 

Statistical bootstrap models [|l], |3|, 11] form particles from clusters of particles, 
and employ the principle of self- similarity. The simplest, "generic", bootstrap 
equation has the form 

p(m) = 5(m — m ) + — / dm 1 ...dm n x 5(m — ■m i )p(mi)...p(m n ), (10) 

where p(m) is the particle spectrum (here, for a moment, mesons and baryons 
are not distinguished). Equation (|T0|) can be nicely solved with help of Laplace 
transforms |3|, |I5|], yielding the asymptotic solution p{m) ~ exp(m/T), with 
T = mo/log(— log -). More complicated bootstrap equations involve integration 
over momenta, more degrees of freedom, different combinatorial factors 0, how- 
ever, irrespectively of these details, they always lead to an exponentially growing 
spectrum. It can be shown, following e.g. the steps of Ref. |16| , that the model 



leads to equal Hagedorn temperatures for mesons and for baryons. This is quite 
obvious. Since baryons are formed by attaching mesons to the "input" baryon, 
the baryon spectrum grows at exactly the same rate as the meson spectrum. 
Specific calculations confirm this simple observation. Thus the bootstrap idea is 
not capable of explaining the different behavior of mesons and baryons in Fig. [I]. 



3.2. Dual String models 

The Dual String models also date back to pre-QCD times. Their greatest 
success is a natural explanation of the Regge trajectories - a basic experimen- 
tal fact which remain a serious problem for other approaches. Similarly to the 
bootstrap models, the Dual String Models lead to exponentially-growing spectra, 
but they do give the demanded effect of T mcson > T baryon , at least at asymptotic 
masses |J. 

Let us analyze mesons first. The particle spectrum is generated by the 
harmonic-oscillator operator describing vibrations of the string, 

oo D 

N = j2T, ka l^ (ii) 
k=ifi=i 
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Fig. 5. Meson and baryon string configurations. 



where k labels the modes and p, labels additional degeneracy, related to the 
number of dimensions [[/J. Eigenvalues of N are composed in order to get the 
square of mass of the meson, according to the Regge formula 

a'm 2 — a = n, (12) 

where a' ~ lGeV -2 is the Regge slope, and «o ~ is the intercept. Here is an 
example: take n — 5. The value 5 can be formed by taking the k = 5 eigenvalue 
of N (this is the leading Regge trajectory, with a maximum angular momentum), 
but we can also obtain the same m 2 by exciting one k = 4 and one k = 1 mode, 
alternatively k = 3 and k = 2 modes, and so on. The number of possibilities 
corresponds to partitioning the number 5 into natural components: 5, 4+1, 3+2, 
3+1 + 1, 2+2+1, 2+1+1 + 1, 1+1+1 + 1+1. Here we have 7 possibilities, but the 
number of partitions grows very fast with n. Partitions with more than one 
component describe the sub-leading Regge trajectories. With D degrees of free- 
dom each component can come in D different species. Let us denote the number 
of partitions in our problem as Pr>{n). For large n the asymptotic formula for 
partitio numerorum leads to the exponential spectrum according to the formula 

p{m) = 2a'mP D {n), P D {n) ~ J^(-£-\ * exp I 2^/~ J . (13) 




2n \24n 

where n = a'm 2 . We can now read-off the mesonic Hagedorn temperature 



^meson — n \ n . ■ (14) 

27T V Doc 

Now the baryons: the "Mercedes-Benz" string configuration for the baryon is 
shown in Fig. [|. The three strings vibrate independently, and the corresponding 
vibration operators, N, add up. Consequently, their eigenvalues rii, n 2 , and n 3 
add up. Thus we simply have a partition problem with 3 times more degrees of 
freedom than in the meson. The replacement D — ► 3D in flHf ) leads immediately 
to 

T baryon = — J — , (15) 
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Fig. 6. Predicions of the scalar string model of Ref. |l0[, with D = 6. 



such that 

-^rneson/ -^baryon \/3. (16) 

We stress that the presented picture is fully consistent with the Regge phe- 
nomenology. The leading Regge trajectory for baryons is generated by the ex- 
citation of a single string, i.e. two out of three numbers n, vanish (this is the 
quark-diquark configuration). The subleading trajectories for baryons come in a 
much larger degeneracy than for mesons, due to more combinatorial possibilities. 
The slopes of the meson and baryon trajectories are universal, and given by a'. 
We stress that the "number-of-strings" mechanism described above is asymptotic. 
Thus, there is a problem in applying string models to the experimentally acces- 
sible range of m. This range is not asymptotic enough to use Eq. (|i~3f ). From the 



Regge formula (12) we find immediately that for m in the range 1 — 2GeV the 
values of n lie between 1 and 4, hence n is not large enough to justify the form 

One can do better by using an improved asymptotic formula, derived in Ref. 
T0| . The results obtained in the scalar string model fl(| are displayed in Fig. 



Here the formula for the meson spectrum is 

Pmes(m) = 36 x p sca i ar (m), p scaiar (m) = — — j^— — —ml u (^—), (17) 

y±KOt Tfll mes ) -'mes 

where I2 is a modified Bessel function, T mes is the meson Hagedorn temperature 
(the only adjustable parameter here), and v = 1 + D/2, with D denoting the 
number of transverse dimensions. The factor of 36 = 6 x 6 is just the spin-flavor 
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q q 
Fig. 7. qqqq and glueball configurations. 

degeneracy of the qq configuration [|l(J . For the baryons we fold the three scalar- 
string densities, p SC aiar(jn)- We use 56 (rather than 36) copies of the string, which 
is the degeneracy of the baryon multiplet in the ground state. We notice good 
agreement with data in Fig. |6], for D = 6. Note that both curves are fitted with 
only one parameter, T mes . For lower values of D one can fit the mesons equally 
well, but too many baryon states are predicted. 



3.3. Exotics as dual strings 



During this workshop we have heard many talks on hadron exotics. If an exotic 
is a multi-string configuration, e.g. as in Fig. then the corresponding spectrum 
will grow exponentially with the Hagedorn temperature inversely proportional to 
the square root of the number of strings. For instance, T q g q g = ^T meson . This is 
reminiscent of the effect described in Ref. |18]. For the glueballs, described by 
the closed string in Fig. (|7|), we get T G = T meson . 

Thus, according to the string model, the qqqq grow more rapidly than non- 
exotic mesons and baryons, and glueballs grow at the same rate as mesons. 



4. Other approaches 

In the remaining part of this talk we will, in a sense, work against our results 
presented in previous sections, where have we argued that the plots of Fig. [I] 
are linear, and offered an explanation of the difference between the mesonic and 
baryonic Hagedorn temperatures within the Dual String Models. 

What if the experimental plots of Fig. [I] are not really linear, and the effect 
of bending down of the curves at higher masses is physical, rather than due to 
incomplete experiments? Below we will show alternative descriptions which do 
not comply to Eq. ([]]), but nevertheless reproduce the present data at least as 
good as the Hagedorn-like fits. 



4-1. Compound hadrons 
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Fig. 8. Fits of the Compound Hadron Model, Eq. (|T9|). 



In the statistical model of nuclear reactions one uses the compound-nucleus 
model 
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In this model the density of states grows at large excitation 



energies, E* , according to the formula 

p(E) ~ (E*) 



where a is a constant. Formula 



-5/4^ 



(11 



can be derived within the Fermi gas model 
More generally, it can be derived in a model where the single-particle orbits 
are equally spaced. One then considers lplh, 2p2h, 3p3h, etc., excitations and 
counts the number of states at a given excitation energy, E*. Amusingly, this 
leads I^TJ to the partitio numerorum formula (pT3|), but now the number n has the 
interpretation n = E*/AE, with AE denoting the level spacing. 

We now use the following Compound-Hadron-Model formula for the mass 
spectra: 

Ae(m-m )exp(2 7 r v /(=l) 
PH = — " (19) 



(m - mo) + (0.5GeV) 2 



where A is a constant, mo is the ground-state mass, and AE is the average 
level spacing. The constant 0.5GeV in the denominator has been introduced ad 
hoc, similarly as in Eq. (|5|), in order for the formula to make sense at m — > m . 
Asymptotically, the power of m multiplying the exponent is —5/4, as in Eq. (|i~8|). 

The underlying physical picture behind compound hadrons is as follows: 
hadrons are bound objects of constituents (quarks, gluons, pions). The Fock 
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space contains a ground state, and excitations on top of it. In the case of the 
compound nucleus these elementary excitations are lplh, 2p2h, 3p3h, etc. states. 
In the case of hadrons they are formed of qq and gluon excitations, e.g. for 
mesons we have qq, qqg, qqqq, qqgg, etc. We can form the excitation energy 
(hadron mass) by differently composing elementary excitations. This bring us 
to the above- described combinatorial problem pl| . It seems reasonable to take 
zero ground-state energy for mesons, m™ es = 0, since they are excitations on top 
of the vacuum. For baryons we take mj" = 900MeV, which is the mass of the 
nucleon. The quantity AE is treated as a model parameter and is fitted to data. 

The results of the compound-hadron-model fit, Eq. ([19]), are shown in Fig. 
S . The curves are slightly bent down, compared to the Hagedorn-like fits of Figs. 
1,|, which is caused by the square root in the exponent of Eq. ([19|) . But the fits 
are at least as good, or even better when the fit region is extended torn = 2GeV. 
Numerically, the least-square fit for m up to 1.8GeV gives AE mes = lOOMeV 
for mesons, and AE ha,r = 106MeV for baryons. The proximity of these numbers 
shows that the scales for mesons and baryons are similar, as should be the case. 

The obtained values for AE mes mean that the corresponding natm = 1.8GeV 
is around 18 for mesons and 9 for baryons. Such values of n are sufficiently large 
to justify the use of the asymptotic formulas. 



4-2. Combinatorial saturation and the light-flavor-desert hypothesis 



There is a possibility of an interesting effect we wish to point out. It is natural 
to expect that a bound hadronic system has an upper limit for the excitation 
energy. It is helpful to think here of bags of finite depth. Thus, in constructing 
the single-particle Fock space for bound objects we should have a limited number 
of quanta to our disposal. If such a limit is put into the Compound Hadron Model, 
it will result in a maximum number of states that can possibly be formed out of 
light quarks . We can call it the "light- flavor-desert hypothesis": above a certain 
mass there are no more light-flavor resonances. Certainly, this is tangential to 
the conventional wisdom that the Regge trajectories should continue indefinitely. 
Note, however, that infinite Regge trajectories have recently been challenged by 
Brisudova, Burakovsky and Goldman, who claim that they should stop around 
m ~ 2.7GeV. Amusingly, this is consistent with the presently-available data. The 
cumulants if Fig. [I] flatten-out in that region. 



4-3. Quark models 

Many talk in this workshop were devoted to variants of the quark model. Here 
we present the result of counting of states in the model of Refs. |22| , E3J, as made 
by Freund and Rosner [f|. 

When we look at Fig. || we again see good agreement in the predicted and 
experimental number of states. This is not at all surprising, since the quark 
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Fig. 9. Experimental cumulants and the predicions of the quark model of Ref. 23 
as counted in Ref. ||, indicated by dots. 



model is designed to fit the data "state by state" in the low-mass regime. As for 
other approaches, spectra at higher m would be needed to verify the predictions. 



5. Final remarks 

There are many fundamental questions which should be cleared when more 
experimental data on hadron resonances are available: Is the Hagedorn hypothesis 
of exponentially-growing spectra indeed correct, or is the growth weaker at higher 
masses? Do the Regge trajectories continue for ever, or stop? Consequently, is 
there a light-flavor desert above a certain mass? Are there exotic states, if so, 
at what rate do they grow?... Certainly, the spectrum above 2GeV may reveal 
many answers and help us to verify various models and approaches. 
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However, even the presently-available spectrum allows for interesting specu- 
lations. Recall the remarks made here by Leonid Glozman, concerning the parity 
doublets in the N and A spectra above 2GeV [24]. Almost all states in that region 
can be paired, and such a regularity suggests that the data in that region may 
be complete! This, in turn, indicates that the bending down of the cumulants in 
Fig. Q may be a physical, rather than experimental effect. 

Another important aspect, not touched in this talk, are the thermodynamical 
implications of the presence of two distinct Hagedorn temperatures for the phe- 
nomenology of heavy-ion collisions, transition to quark-gluon plasma, etc. This 
will be discussed in f2lfl . 

The author thanks Keith R. Dienes for many profitable e-mail discussions on 
the issues of hadron spectra in string models, as well as to Andrzej Bialas, Andrzej 
Horzela, Jan Kwiecihski, and Kacper Zalewski for numerous useful comments and 
encouragement . 
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